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We investigate the upper critical field in a dirty two-band superconductor within quasiclassical Us- 
adel equations. The regime of very high anisotropy in the quasi-2D band, relevant for MgB2, is con- 
sidered. We show that strong disparities in pairing interactions and diffusion constant anisotropies 
for two bands influence the in-plane H C 2 in a different way at high and low temperatures. This causes 
temperature-dependent H C 2 anisotropy, in accordance with recent experimental data in MgB2. The 
three-dimensional band most strongly influences the in-plane H C 2 near T c , in the Ginzburg-Landau 
(GL) region. However, due to a very large difference between the c-axis coherence lengths in the 
two bands, the GL theory is applicable only in an extremely narrow temperature range near T c . 
The angular dependence of H C 2 deviates from a simple effective-mass law even near T c . 
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I. INTRODUCTION 

There is a strong evidence of the multigap nature of 
superconducting state in the recently discovered^ com- 
pound MgB2. The concept of multiband superconduc- 
tivity was introduced in 2 * 3 - for the case of large dis- 
parity of the electron-phonon interaction for the differ- 
ent Fermi-surface sheets. For MgB2, first-principles cal- 
culations of the electronic structure and the electron- 
phonon interaction2«ii£iL2ai have revealed two distinct 
groups of bands, namely strongly superconducting quasi- 
two-dimensional er-bands and weakly superconducting 
three-dimensional 7r-bands. Quantitative predictions 
for various thermodynamic and transport properties of 
MgB2 were made in the framework of the two-band 
model, 10 . 11 ! 12 . 13 

A large number of experimental data, in particular 
tunnclingfi^ii^ point contact measurementsjiSiiLiS and 
heat capacity measurements^ directly support the con- 
cept of a double gap MgB2. Intraband impurity scatter- 
ing in both bands may vary in large limits, while intcr- 
band scattering is always weak due to the disparity of 
a- and 7r-band wave functions^ This explains the ex- 
tremely weak suppression of T c by impurities and the 
weak correlation between T c and the resistivity. There- 
fore, a unique feature of the MgB 2 is that the two-gap na- 
ture of superconductivity persists even in the dirty limit 
for the intraband scattering rates. 

Superconductivity in the two bands is characterized 
by different energy and length scales which show up in 
several properties of a superconductor. Particularly in- 
teresting are the properties of the mixed state. The c- 
axis Abrikosov vortex structure in MgB2 was studied by 
STM in Ref.Hl which probes mainly the weakly super- 
conducting 7r-band. A large vortex core size compared to 
estimates based on H C 2 and the rapid suppression of the 
apparent tunneling gap by small magnetic fields has been 



reported. These observations can be naturally explained 
within the two-band modeliSL 2 ^ 

One of the most spectacular consequences of the 
two-band superconductivity is the unusual behavior of 
anisotropy factors for different physical parameters^ 3 - 
ft was demonstrated that in clean MgB2 samples the 



anisotropy of the London penetration depth 



24.25 



7a 



has 



to be very different from the anisotropy of the upper criti- 
cal fieldrL^l -y c2 . Both anisotropy factors should strongly 
depend on temperature and have opposite temperature 
dependencies: 7a is expected to increase and 7 C 2 is ex- 
pected to decrease with temperature. Strong tempera- 
ture dependence of 7 C 2 has been reliably confirmed by 
experimenti 2 ?. 2 ?' 30 . 31 . 32 . 33 Typically, 7 C 2 drops from 5-6 
at low temperatures down to ~ 2 near T c . 

In this paper we consider in detail the behavior of the 
upper critical field for different field orientations for the 
case of a dirty two-band superconductor with weak inter- 
band scattering. The model is based on the multiband 
generalization of the quasiclassical Usadel equations. 34 
The same model has been used recently to describe vor- 
tex core structure in MgB 2 . ■ — The general equations for 
determination of the upper critical field within this model 
have been derived in recent paper l3ol However, calcula- 
tions in this paper have been done only for the case of 
small band anisotropies. In this paper we address the 
case of very high anisotropy in the quasi-2D band, more 
suitable for MgB2. 

We demonstrate that the strong temperature depen- 
dence of the ir c 2-anisotropy exists also in the dirty case 
and therefore represents a general property of a two-band 
superconductor. The main reason for this dependence is 
the strong reduction of the in-plane upper critical field by 
the weak 7r-band in the very narrow temperature region 
near T c . This also leads to the significant upward curva- 
ture of the temperature dependence of the in-plane upper 
critical field near T c . This behavior illustrates breakdown 



2 



of the anisotropic Ginzburg-Landau (GL) theory for de- 
scription of this superconductor. We demonstrate that, 
due to the large difference between microscopic coher- 
ence lengths in the c-direction for the two bands, the 
anisotropic GL theory is applicable only within the ex- 
tremely narrow temperature range near T c . 

We analyze the angular dependence of the upper crit- 
ical field and show that it strongly deviates from the 
standard "effective-mass" dependence predicted by the 
anisotropic GL theory. Contrary to naive expectations, 
these deviations are strongest for temperatures quite 
close to T c (at T ~ 0.9T C ) and vanish only for temper- 
atures extremely close to T c (for (T c - T)/T c < 1%). 
In the past the angular dependence of the upper critical 
field have been studied in Ref. 40 for a clean two-band 
superconductor. It was shown that for the case of two 
weakly deformed spherical Fermi surfaces with opposite 
anisotropies the angular dependence also strongly devi- 
ates from the "effective-mass" law. 

The paper is organized as follows. In section [H] we 
present Usadel equations for a two-band superconduc- 
tor and introduce parameters relevant for MgB2. In sec- 
tion [ffi] we derive equation for the upper critical field in 
the c-direction and obtain the exact asymptotics at small 
and high temperatures. In section Hvl we consider the in- 
plane upper critical field. We derive general equations for 
determination of this field and study solutions of these 
equations in different regimes. We demonstrate that the 
GL result for the in-plane H c2 is valid only within a very 
narrow range of temperatures. We also numerically cal- 
culate in-plane H c2 and the anisotropy parameter "f c2 in 
the whole temperature range. In section we study the 
angular dependence of the upper critical field and ana- 
lyze quantitatively the deviations from the effective-mass 
law. 



II. THE MODEL: USADEL EQUATIONS FOR A 
TWO-BAND SUPERCONDUCTOR 

We consider a two-band superconductor with weak 
interband impurity scattering and rather strong intra- 
band scattering rates exceeding the corresponding energy 
gaps (dirty limit). In this case the quasiclassical Usadel 
equations^ are applicable within each band. The mixed 
state in this case is described by the system of coupled 
Usadel equation s 21 ^ 4 



T> a , 



3 

— <A a G a , 

0,n 



G a (V 3 ~^-A 3 ) 2 F a -F a Vp a 



(la) 
(lb) 



where a = 1,2 is the band index, j = x,y,z is the co- 
ordinate index, A is the matrix of effective coupling con- 
stants, T> a j are diffusion constants, which determine the 



coherence lengths £ a j = y/T> a j /2nT c , G a , F a and A a 
are normal and anomalous Green's functions and the pair 
potential, respectively, and ui = 2irT(s + 1/2) are Mat- 
subara frequencies. Bearing in mind the application to 
MgB2, in our notations index 1 corresponds to cr-bands 
and index 2 to 7r-bands. All bands arc isotropic in the xy 
plane, T> ax — T> ay and anisotropic in the xz plane with 
the anisotropy ratios j a = ^/T> ax /T> az . The multigap 
Usadel equations for general case, taking into account 
also interband scattering, have been recently derived in 
Ref.lM 

The selfconsistency equation can be rewritten in the 
form 



W^i- W 12 A 2 =2ttTY / (fi- — ) + Ailn^p, (2a) 



~IF 2 iA 1 + VF 2 A 2 =2^T^(f 2 - ^ij + A 2 lnp, (2b) 
with the following matrix W a p 



„, -A + \M 2 + A 12 A 2 i A + y/A 2 + A 12 A 2 i 

Wi = — , W 2 = — 



Det 

W 12 = Au/Det, TU 2 i = A 2 i/Det, 



Det 



(3) 



A = (An - A 22 )/2, Det = A n A 22 - A 12 A 2 i, W{W 2 = 
W 12 W 21 . 

The electron-phonon interaction in MgB 2 was calcu- 
lated from first principles in a number of papers&iSiii 
Here we use the effective coupling constants A^ from Ref. 
[U An » 0.81, A 22 w 0.278, A i2 » 0.115, A 2i w 0.091, 
from which we obtain values of W a p used in numerical 
calculations, 

Wi« 0.088, TU 2 w2.56, W 12 ^ 0.535, W 21 ^ 0.424. (4) 

The relative role of the weak band is characterized by the 
ratio S\ 2 = W\/W 2 ^ which in the case of MgB 2 is rather 
small, Si 2 ~ 0.034. This ratio will be used below as a 
small parameter in our model to derive various approx- 
imations for the upper critical field. Another important 
small parameter is the ratio of diffusion coefficients in the 
(T-band, T>\ z jT>\ x . We will show in this paper that these 
two parameters, «Si 2 and T>i z /T>i Xl influence differently 
H c2 for parallel field at high and low temperatures thus 
causing the temperature dependence of the anisotropy. 

In the following we consider separately the cases when 
the field is parallel and perpendicular to the ab-plane. 



III. FIELD IN THE c-DIRECTION 

Let us first study the case when the magnetic field is 
oriented along c-axis. The upper critical field is deter- 
mined by the linearized Usadel equation 



ivF a 



-V 2 F 

v x r a 



2-kHx 



(5) 
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and selfconsistency equations Solving these equa- 

tions, we arrive at the equation for (symbol _L de- 
notes the field direction perpendicular to the (ab)- plane) 



The well-known asymptotic solutions of this equation at 
low and high temperatures are respectively 



hi- 



w 2 



y \ tH 2 



(6) 



where t = T/T c , H a = 2T C ^„/V ax , g(x) = 0(1/2 + x) - 
■0(1/2), and ip(x) is a digamma function. We also obtain 
a relation between Am and An? near H, 



L r2 



A 2 = 



W 2 iA 01 



W-2 



14 



■9 



(ml) 



(7) 



In the absence of coupling to the weak 7r-band (W% = 
0) or in the case of identical diffusion constants (T>i x = 
T>2 X ), the upper critical field H s c2 is given by the standard 
Maki - de Gennes equation 38-39 



Hl/t)=g[H s c2 /(tH 1 )}. 



(8) 



e"^/ 4 - 0.140, t < 1, 
2(1 - <)/tt 2 w 0.203(1 - t), 1 - t < 1, 



(9) 

where ~ 0.577 is Euler constant. In the temper- 
ature range near T c one can obtain from Eq. JBJ) the 
following simple expression for for arbitrary ratio 
S 12 = W 1 /W 2 : 



H^2 = 2(l + S 12 )(l-t) 
7T 2 (1 + S 12 V 2x /V lx 



(10) 



At small temperatures, T <C T c , Eq. 10 also has an 
exact solution (see also Ref. I35T) 



H c2 (0) -i/ c s 2 (0)exp 



W^ + W2-In(rx) , /(Wi + W 2 -ln(r x )) 2 



Wiln (r x ) 



(11) 



with r x = T>i x /T> 2x . For MgB2 the parameter W\ 
is small and typically the inequality W\\a.{r x ) <C 
(W 2 — In (r x )) 2 /4 is valid. In this case we can expand 
Eq. [fTTJ) with respect to W\ and obtain a much simpler 
result 



FF 2 - In {V lx /V 2x ) 



The 7r-band strongly influences the upper critical field 
only if it is very dirty, T> 2x <C Din exp(— W 2 ). In this 
limit we obtain 3 ^ 



at T c and H c2 (0): 



dH, 



c2 



dH. 



c2 



dT 



dT 



l + S 



12- 



",•>(» 1 = //;,(()) j 1 + 5ia In ^ 



(13a) 
(13b) 



The signs of the above corrections to the universal curve 
following from Eq. (JHJ are positive if D 2x < D\ x and 
negative for D 2x > D\ x . 



IV. FIELD IN THE a -DIRECTION 



# c2 (0)«i4 2) (0)exp (-W 2 ) 

with H%\0) ee (exp(-7E)/4)£r 2 . 

For the case W\ <C W 2 realized in MgB2, the upper 
critical field is typically determined by the strong band 
(except for the limit of very small diffusivity T> 2x in the 
second band). A small correction due to the weak band 
can be found from Eq. (jSJ using an expansion with re- 
spect to the small parameter £12 = W\/W 2 . In particu- 
lar, we found very simple expressions for the slope of H c2 



A. General relations 

The upper critical field in the a-direction (= y- 
direction) is determined by the linear equations for the 
Green's functions F a in two bands 

uF a - — V x F a + — I —j— I F a = A a (14) 

with u> = 2ttT(s+1/2) and the self-consistency conditions 
(J2J). A technical difficulty of this problem is that, due to 
the difference in the anisotropy factors for the two bands, 
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7 Q , the harmonic oscillator operators in Eq. I|14f) have 
unmatching sets of eigenstates. We will use an expansion 
with respect to the eigenfunctions (Landau levels) of the 
strong [1st] band, ^> n (x), which are defined as solutions 
of the oscillator equation 



(15) 



In particular, the eigenvalues e n and ground state eigen- 
function are given by 



En = v^i^teS- (n + 1/2) , 



*o(a:) 



2H 

7i $o 



1/4 



exp 



irHx 2 
7i $o 



(16) 
(17) 



where 7 Q = \jT) ax /T> az are the band anisotropies. In 
the case of MgB2 the first band is quasi- two-dimensional, 
i.e., 71 ^> 1,72- Substituting expansions 



A Q (a;) =^A a , n * n (a;); F a = ^F a ^ n {a 

n n 

into Eq. Ill4|l . we obtain 

oo 



^-^2,71 + e,im '' •' 

m— 



CJ + 

A 2 .„ 



(18a) 
(18b) 



with e. 



22, ( 2nH\ 

2 \ *o j 



„2 ©2 



2 VI 



The only 



nonzero matrix elements e„m are at m = n and to = n±2: 



7Ti7 , / 1 

yV 2x T> 2z [n + - 



72 7i 
7i 72 



(19a) 



^n — 2,n — €n,n— 2 



irHs/V 2 ~^D 2 ~ z \Jn{n - 1) / 72 71 



71 72 



(19b) 



Neglecting the small ratio 71/72 in comparison with 
72/71 we obtain 



n + - l u?2, 



^n — 2.n ^n,n — 2 



^n(n - 1) 



-io 2 



(20a) 
(20b) 



with u>2 = (~kH /^ofD 2 Z 7i. This approximation for the 
matrix elements is equivalent to the local approximation 
for the F-function in the 7r-band described in Appendix 
1X1 Therefore, we can rewrite the equation for F 2n as 



W-F2n+ e n,n-2-F2,n-2+en,n-p2,n+en,n+2-p2,n — A2, n . (21) 



At n = the term e n ,n—2F 2in — 2 has to be skipped. This 
means that even Landau levels, n = 2i, do not mix with 
the odd Landau level, n — 2i + 1. For calculation of the 
upper critical field it is sufficient to consider only even 
Landau levels. The self-consistency equations in terms 
of the expansion coefficients are given by 



WiAi, n - W 12 A 2 , n = 2irT ^ ( — ] A Un + A hn In 1 



UJ>0 



UJ + E n U 



-W 21 A hn + W 2 A 2 , n = 2nT J2 ( F 2,n - —) + A 2 ,„ In i. 



w>0 



(22a) 
(22b) 



To simplify further analysis we introduce the reduced variables 

z = uj/w 2 = t 2 (s + 1/2), F 2>i = w 2 F 2j2i 
with t 2 = 2irT/w 2 = 2&qT / (HT) 2z ji) and s is the Matsubara index. Then equations for F 2i i(z) are given by 

z+^\ F 2 , -(l/V2)F 2il = A 2fi , 



-y/i{i- l/2)F 2ii _i + (z + 2i + 1/2) F 2 , 4 - y/(i + l)(z + l/2)F 2 ,i+i = A 2 , 2 *. 
The formal solution of Eq. I|23|l is given by 

oo 

h,i(z) =^A M (z)A 2 , 2 ,, 

J=0 



(23a) 
(23b) 
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where the matrix Ai (z) is defined as solution of equations 



-V*(*-l/2)Ai_u + (z + 2 2 + 1/2) i4ij - ^/(i+l)(i+l/2)A+i,i = 

I 



(24a) 
(24b) 



Using this solution we represent the self-consistency con- 
ditions for even Landau levels in the form 



WiAi,ai-Wi 2 A 2 , 2< =Ai,2i [in- -g ' /7(4 ' " 



t# 1 



(25a) 



v - 

-W A 2lAi, 2i +W r 2A 2 ,2i=X] C/ i > j(*2)A2,2i+A 2)2i In - 



J=0 



with 



2T e $o 



(25b) 



(26) 



where symbol || denotes the field direction parallel to the 
{ah)- plane, and 



^ij(*2)=*aX] ( [ta(* + l/2)]- 



t 2 (s + l/2) 



(27) 



We again used notations t = T/T c and = 

if) (1/2 + x) - -0(1/2)- We show in Appendix |X| that 

Uij(t 2 ) can also be related with the oscillator matrix 
element of the function g [x 2 /t 2 ] 

/oo 
dxip 2i (x)i/; 2 j(x)g [x 2 /t 2 ] , 
-OO 

exp (-X 2 /2) H n (x) 

Wn(x) — " ; ; , 

where H n (x) are Hermite polynomials. 



B. Temperatures not close to T c . High-field 
approximation in the 7r-band 

The overall behavior is determined by the value of di- 
mensionless parameter t 2 , which depends on field and 
temperature. To evaluate this parameter we represent it 
in the form 



V 2z H 



(28) 



Because T>i z <C T> 2z and at low temperatures H < H { , 
the parameter t 2 is much smaller than unity almost in 
the whole temperature range except a very narrow re- 
gion near T c . The parameter t 2 becomes of the order 
of one only at (T c - T)/T c ~ V lz /V 2z < 1. Outside 
this region one can replace summation with respect to 
the Matsubara index s in Eq. I|27|l by integration, which 
allows us to reduce it to the following form 



Ui,j(t 2 ) m f itj + (lnt 2 - j E - 2 In 2) S h3l 



where 



dz A, j (z) 



z + 1 



(29) 



= -4/ dxip 2l (x)i) 2j (x)\n(x) (30) 
Jo 

is the universal matrix of constants (in particular, /o^ = 
je + 2 In 2 f» 1.96). Using this representation, we trans- 
form Eq. i|25b(l to the form 



-W 21 A li2l + W 2 A 2 , 2i 



3=0 



1 :,; 1 irk. 



•7 B + 21n2 A 2 , 2l 



(31) 



We refer to this approximation as the high-field regime in the 7r-band. The last equation in combination with 
Eq. Q25al ) determines the upper critical field along a-direction within the "high-field in the 7r-band" regime, at 
(T c — T)/T c ^S> T>i z /T> 2z . Note that in this approximation the temperature dependence exists only in Eq. l|25afl . 
Therefore, once computed, matrix fej allows us to calculate the temperature dependence of H^ 2 in a wide temperature 
range. 

Excluding Ai 2 i 



±l,2i 



Wis 



Wi 



H(ii+1) 



A 



2,2i 



(32) 
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we also derive equations containing only A 2j2 ^ 



-W 2 In i 



V 



g(4i+l) 

, H V 2z \ " * V ~* 4 y \ tH. 

M^J+7, + 21n2 



In I - o ( g ( 4 ' +1 ) 



r 



A 2 ,2i-5Z/i.jA 2 ,2j = 0. 



(33) 



The upper critical field _ff = iJ^, 2 is given by the maxi- 
mum root of the determinant of this linear system. An 
approximate solution can be obtained neglecting coupling 
to the higher Landau levels in the self-consistency equa- 
tions leading to the following equation for H^ 2 



In- 



H 



c2 



tHl 



(34) 



Since Wi <C 1, the right hand side of Ea. l(3^jl is small. As 
a result, in the limit of small t 2 the parallel critical field 



is close to the solution of the Maki - de Gennes equation 
(JSJ) with the effective parameter Hi replaced by H\ from 
Eq. (|26|) . A small correction from the weak band can be 
estimated at low temperatures 



Wi(ln (ZWPis) - 1.96) 
W 2 +ln {V 2z /V lz ) - 1.96 



^(0)^ £ 1(0) i 



(35) 

withff c s J(0H (exp(- 7£ )/4)i/f. 

Combining Eqs. i|12|) and (|35|l we obtain an estimate 
for the anisotropy factor 7 c2 (T) = H l } 2 (T)/ H^(T) at low 
temperatures 



7c 2 (0)«7i 1 + 



W x ln (V 2x /V lx ) Wi (In (V 2z /V lz ) - 1.96) 



^2 + In (V 2x /V lx ) W 2 + ln {V 2z /V lz ) - 1.96, 
As follows from this equation, the anisotropy of H c2 at T = is very close to the anisotropy of the first band 



7c2 (0) W 71 EE y/V lx /V U 



(36) 



(37) 



To estimate the ratio T>x x /T>x z — v f1x ti z /v f1z ti x for MgB 2 , we take v Flx /v Flz ~ 40 provided in Ref . fl2l and assume 
isotropic scattering t\ x rs t\ z . This gives \jT>\ x jT>\ z rs -\/40 » 6.3, which is consistent with the experimental data on 
the H c2 anisotropy in MgB 2 single crystalsi 28 i 29 i 30 i 31 i 32 

C. Ginzburg-Landau region 

In the close vicinity of T c (exact criterion will be established below) one can solve Eq. I|14fl using the gradient 
expansion 

Substituting this expansion into the self-consistency conditions and using relation 27rT^ w>0 (l/cj 2 ) = tt/AT, we 
obtain 



WiAi - W 12 A 2 



{2nHxY \ . 1 
A,l - A, In - 

-Wax Aj + W 2 A 2 = - f -CLV 2 A 2 + ( ^lE^] a 2 ] + A 2 In 1 



-^V^A 1+ ef 2 J A, ) - A, h, 

/2tT#£\ 



(38a) 
(38b) 



with £^ = irT> a i / (8T) . Near T c we can look for solution where S 2 is a small correction, for which we obtain from 
for A 2 in the form 



Woi 

A2 -lf A ^ 2 < 



Eq. J38b 



-Wi 2 S 2 



(1-kHxX 



-£ e V£A 1 +&(- i — 1 Ail+Axln-. 
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Substituting this result into Eq. I|38a|) we obtain the lin- 
ear Ginzburg-Landau (GL) equation for Ai 



i 



A 1 -A 1 ln-=0, (39) 



in which the averaged coherence lengths, £i with i = x,z, 
are defined as 



1 + S 



12 



(40) 



From this equation we immediately obtain the usual GL 
result for the upper critical field at T — ► T c 



ff C 2 = 



$o(l-t) 

' 



(41) 



For comparison with numerical results at lower temper- 



ature we also provide Hjg in units of ii{ 



_ 2y/V lx V lz (1 + Si 2 ) (1 - t) 

iff n 2 y/(T>i x + S 12 V 2x ) (V lz + S 12 V 2z ) 

:{l-t) 



(42) 



ttV 1 + Sx2V 2z /T> x , 

for Wi < W2 and D lx ~ £> 2:c . 

Due to the strong inequality T> 2z T>\ z , in the vicin- 
ity of T c the three-dimensional band strongly reduces the 
upper critical field. This reduction leads to a strong tem- 
perature dependence of the H c2 anisotropy, 7 c2 - 

Let us compare anisotropy parameters at low T and 
near T c . According to Eq. (32J, the anisotropy of H c2 at 
low temperatures is close to the anisotropy of the c-band, 
7c2(0) ~ y/T>i x /T>i z , while the anisotropy ratio near T c 
follows from Eqs. fTOjl and l|l2"} 



lc2(T c ) = JGL — 714/ , (43) 

y J- + o\ 2 u 2z i U\ z 

_ 7i 

Thus the ratio 7 C 2(0)/7c2(^c) is roughly given by 
7c 2 (0) 



lc2{T c 



y/i + Si 2 V 2z /V lz 



(44) 



The larger is the ratio of transport constants, T>i z /T> 2zi 
the stronger is the suppression of "f c2 (T) with increasing 
temperature. 

We obtain now the applicability criterion for the GL 
expansion. Typical scales of the order parameter vari- 
ation near T c are given by the GL coherence lengths 
£,f L {T) = &/V1 - t, with i = x,y and & given by Eq. 
I|4()jl. The GL expansion is valid until the GL coherence 
lengths are larger than the corresponding microscopic co- 
herence lengths in both bands, £,f L (T) > £ a) j. Because 



of the strong inequality £i. 2 <C £2,2, the most sensitive 
condition is 

<g L <X)>b„ (45) 

leading to the following condition for the GL temperature 
range 



T c -T /& 

< max I 72~, A12 

\S2z 



(46) 



Because <C £, 2z and S' 12 1, the applicability of the 
GL approach is limited to an extremely narrow temper- 
ature range near T c , i.e., the situation is very different 
from usual single-band superconductors. The compari- 
son of the GL asymptotic and with the exact solution is 
shown in Fig. ^ where the narrowness of the GL region 
is demonstrated in the inset. 



D. Numerical solution in the whole temperature 
range. 

In the whole temperature range, for an arbitrary value 
of the parameter t 2l the problem can be solved numeri- 
cally. The solution consists of three steps: (i) the matrix 
Aij(z) has to be found from Eqs. I|24|) for the series of 
reduced Matsubara frequencies z — t 2 (s + 1/2), (ii) the 
matrix Uij has to be computed by summation over Mat- 
subara indices s l|27|l and (iii)the upper critical field has 
to be found as the maximum root of the determinant of 
the linear system represented by Eqs. I|25a|) and l|25b(l . 
Due to fast decrease of the nondiagonal matrix elements 
Uij for — 3> 1 , sufficient accuracy is achieved for 
dimension of the matrix less than 30. The result of cal- 
culation of the parallel upper critical field is shown in 
Fig. where the ratio T> 2z jT)\ z — 100 relevant to MgB 2 
was used. Note that when plotted in reduced units, the 
deviations of both ratios n\ 2 l 'iff and H^/H^~ from the 
universal single band curve are small (except from the 
region near T c , in the GL region), in accordance with the 
above discussion. However, one should keep in mind the 
large difference in magnitudes of the characteristic scales 
iff and H^. 

Numerically calculated temperature dependence of the 
anisotropy factor for several ratios T> 2z /T>i z is shown in 
Fig. The anisotropy ratio drops with the increase of 
temperature, in accordance with the estimate l|44l) . This 
result agrees qualitatively with recent measurements of 
temperature-dependent anisotropy in MgB 2 . 28 - 29-30 ' 31 ' 32 
In experiment the change in anisotropy typically is dis- 
tributed over wider temperature range than it is sug- 
gested by the theory. 



V. TILTED FIELDS 

The upper critical field for magnetic field tilted at angle 
with respect to z axis in [zy) plane is determined by 
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FIG. 1: Temperature dependence of the upper critical field in 
the a-direction normalized to h\ defined in Eq. 1261 . Single- 
band curve is normalized to the corresponding scale H\ . Solid 
circles show the dependence obtained within the "high-field in 
the 7r-band" approximation (Eqs. 129H and 13ft ). Inset: com- 
parison between the exact solution and the GL asymptotics 
Eq. JH. 



C\J 

^° 5 



2z 1z 

— 600 

— 300 

—□—60 




0.0 



0.2 



0.4 0.6 
t 



0.8 



1.0 



FIG. 2: Temperature dependence of the anisotropy of the 
upper critical field for several ratios T>2z/'D\ Z . 



the coupled linear equations for the Green's functions F a 
in two bands 



ujF n 



with 



2 I $o j 



F a = A a (47) 



V a (0) = V ax cos 2 + V az sin 2 9 



(48) 



and the self-consistency conditions @. 

Therefore the 7f c2 -problcm of the upper critical field in 
tilted field reduces to the in-plane iJ C 2-problem by sub- 
stitution T> az — > T> a (0). It is convenient to introduce the 
angular-dependent anisotropy parameters 



l a {9) 



V a {6) 



7 2 cos 2 1 



(49) 



sin 



Such defined anisotropy parameters vary from 1 to 7 Q 
when angle varies from to tt/2. 

Following the route of the previous Section, we again 
use expansion with respect to the Landau levels of the 
strong band, defined by Eq. JT3 with V lz -> T> 1 (9). 
The ^-function of the strong band is given by 



Fin — 



e n (0) 



with the eigenvalue 

Sn(0) -- 

Ht{9) -- 



2vrT 



H 



m(0) 

2T e $ 



(2n + l), 



The matrix elements for the harmonic oscillator operator 
of the weak band are given by 



e nn = —V 2 (0) ll (0) (n 



72 (ft) 

71 (ft) 



2irT , ^ f 1 



£n,n-1 = ^ " 7, V 2 (9)l\ (ft) 1 ~ 



$0 



7i (0) 



W) { 7()) 2 



with 



h(9) 

a 7 (0) 



2T$ 



2T$o\/ cos26 ' + 7] "sin 



HV 2 {9)n(9) HV 2x (cos 2 < 
72(9)\ 2 _ l + 7r 2 tan 2 6» 



7~ 2 sin 2 0) 



li(0) 



1 + 72 tan 2 



9 



Note that at arbitrary angle we can not use inequality 7 2 (#)/7i(#) "C 1 any more. The system of equations for the 
reduced ^-function at even Landau levels, F 2n = (2TrT/t 2 (9))F 22i , a t arbitrary tilt angle is given by 

-(1 - a 7 ) - l/2)A,<-i + (z + (1 + « 7 ) hi F 2 ,i - (1 - a 7 )^(i + l/2)(i + l)F 2 , i+1 = A 2 , 2i (50) 



with z = t 2 (9)(s + 1/2). 

At small tilt angles, 9 <C 1, one can solve Eq. lj5U|l using 
perturbation theory with respect to 2 . The quadratic 
angular correction can be obtained neglecting coupling to 
the higher Landau level. This leads to equation similar 
to Eq. 10 with replacements 



Hi 



Hi 



Hi(0) 
H 2 {9) 



2T C $ 



V 2 (9h 1 (6)(l + o h y 



At small angles we obtain quadratic in 9 corrections to 
typical fields 



i + (i- 7 r 2 )- 



H 2 (9)^H 2 l + (l-7 2 ~ 2 )- 



At low temperature one can derive an exact formula for 
small-angle correction 



H c2 (9) - H c2 (0) Kj PL_l 1 7f a +7 -2 _ (7 2 " 2 -7r 2 )(^ 2 -^ 1 -lnr 3; ; ^ 
Hc2(0) 2 V 2 V 1 2 \/(W 2 + W 1 -lnr x ) 2 + 4W 1 \nr x 

with r x = T>i x /V 2x . In the case of small correction from the weak band, 4Wi \nr x <C (W2 — lnr x ) 2 , we obtain a 
simpler formula for 9 <C 1 

fl*(6) 71 + (^-Inr.) 2 J' (52) 

For parameters of MgB 2 this formula gives an estimate almost identical to the exact result. 

At large tilt angles, cos# <C 1, inequality r ) 2 {9) -C 7 i(0) is restored and we can utilize the approximations used 
for the case of in-plane field. In particular, at low temperatures the approximate angular dependence is given by a 
formula similar to Eq. (|35|l . 

g g(o) / M'-^w-'-^ m 

/cos2 + 7r 2 sin 2 I ^+ln PljCOta a ' 9+Plz -1.96 / 



r 



In the whole angular range we calculated the upper 
critical field numerically following the procedure outlined 
in Sec. lIVDl As input parameters we have used the val- 
ues 71 = 6.325, 71 = 0.816 which follow from the elec- 
tronic band-structure calculations in MgB 2 . We have 
also used the relation T>\ x = 0.2T> 2x - the reason for this 
choice was discussed in Ref. El The examples of the cal- 
culated angular dependence for T/T c — 0.1 and 0.95 are 
shown in Fig. [3] . We also show fits to a simple effective- 
mass law, routinely used to describe angular depen- 
dence of H c2 in anisotropic superconductors, H c2 {9) = 

H c2 c / \ /cos 2 9 + j~ 2 sin 2 9. Due to the contribution from 



the 7r-band, one can see significant deviations from this 
law at high temperature. To enhance these deviations we 
plot in Fig. 0] the angular dependence of the combination 
A{9) = {H c2 , z {9)/H c2 , c f + {H c2 , x {9)/H c2 , a f for several 
temperatures, (for the effective-mass law A{9) = 1 for 
all 9). We find that always A{9) < 1 and the maxi- 
mum deviation from unity is achieved around 9 ~ 74°. 
At high temperatures one can derive a very simple for- 
mula for A(9) at small angles, 9 <C 1, A(9) f» 1 — 
(1/7gl — l/7c2)^ 2 - Quantitatively, the deviations from 
the effective-mass law can be characterized by the param- 
eter (L4 ma x = maxg |1 — A(6)\. Fig.[5]shows the tempera- 
ture dependence of this parameter. At low temperatures 
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FIG. 3: Examples of angular dependence of the upper critical 
field at low and high temperatures. Fits to the effective-mass 
dependence are also shown. 



FIG. 5: Temperature dependence of the parameter (5*4 max = 
maxg 1 1 — A{9) | characterizing deviations from the effective- 
mass dependence of H C 2- Inset shows the dependence of this 
parameter on the i/ C 2-anisotropy. 
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FIG. 4: Plots of the parameter A(6) = {H c2 , z {8)/H c2 , c ) 2 + 
(H C 2,x(9) / H C 2,a) 2 vs sin 2 8 at different temperatures revealing 
deviations from the simple effective-mass law. Left panel: 
temperatures not very close to T c , right panel: temperature 
region near T c . 



deviations from the effective-mass law are at the level 
of several percents. These deviations progressively grow 
with the temperature reaching ~ 19% at T/T c = 0.92 and 
then rapidly decrease when the temperature approaches 
the narrow GL region near T c . As the deviations from the 
effective-mass dependence have exactly the same origin 
as the temperature dependence of the anisotropy, it is in- 
teresting to correlate these deviations with the anisotropy 
change. Inset in the Fig. [5] shows plot of the parameter 
<5-4max vs the i? C 2-anisotropy. One can see that the dis- 
tortion of the angular dependence is maximum when the 
anisotropy is approximately at the midpoint between the 
low-tempcrature and GL limits. Experimentally, it was 
found that the angular dependence of H C 2 in MgB2 in- 
deed deviates from the effective mass law^ii and the 
shape of these deviations qualitatively agrees with our 
calculations. 



VI. CONCLUSIONS 

We have calculated the upper critical field in a dirty 
two-band superconductor within the quasiclassical Us- 
adel equations, bearing in mind the regime of very high 
anisotropy in the quasi-2D band relevant for MgB2. Fol- 
lowing Ref. 

El 

we have assumed that the interband scat- 
tering is negligible even in the dirty limit in both bands. 
Most of MgB2 samples are in dirty limit, except for single 
crystals, where the dirty limit conditions are fulfilled in 
the 7r-band but not fulfilled in the er-band. 41 Still, as ar- 
gued in Ref. 21, our results should be qualitatively appli- 
cable to MgB2 single crystals, if one considers the coher- 
ence length £i as a phenomenological parameter instead 
of expressing it via the diffusion constant T>±. 

We have considered separately the cases when the 
field is parallel and perpendicular to the basal plane. 
We have found that at low temperatures both critical 
fields are mainly determined by the strong band and only 
weakly deviate from the universal Maki - de Gennes re- 
sult. The low temperature anisotropy is mainly deter- 
mined by the anisotropy of diffusion constants in a quasi- 
two-dimensional band. However, the anisotropy is sup- 
pressed at high temperatures. The reason is that there 
are two important parameters, anisotropy of pairing in- 
teraction and of diffusion constants, which enter the ex- 
pression for the parallel H C 2 in a different way at high 
and low temperatures. This property can be expressed 
as the anisotropy of coherence length £ a b/£c which de- 
creases with increasing temperature. This effect is in ac- 
cordance with the experimental data in MgB2. Note that 
the anisotropy of the penetration depth X c /X a b increases 
with increasing temperature, which is another man- 
ifestation of the two-band model. 

We have also studied quantitatively the dependence 
of case of H C 2 on the angle between the a6-plane and 
the magnetic field direction. Approximate relations for 
H C 2 dependence on titled angle are derived for small and 
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large angles. In the whole angular range numerical calcu- 
lations are performed. The results demonstrate the de- 
viation from the effective-mass dependence. This means 
the breakdown of anisotropic GL theory. Further, we 
have shown that the temperature range of applicability 
of the GL theory is extremely narrow in the considered 
two-band case. 

Another issue is strong coupling corrections to H C 2- 
In this paper the weak coupling approach was used. 
On the other hand, it is known from work on 
isotropic superconductors^ that strong coupling correc- 
tions renormalize the absolute value of H C 2 by the factor 
(1 + A) Q , where A is the coupling constant and a ~ 2. 
Since electron-phonon coupling in MgB2 is relatively 
strong (according to Ref. I 111 An ~ 1), these corrections 
are important for calculation of absolute values of H C 2- 
However, we do not expect qualitative changes in the 
temperature and angle dependencies of the anisotropy 
ratio calculated in the present paper. Extension of our 
results to the strong coupling Eliashbcrg regime is an 
interesting subject for future work. 

We acknowledge valuable discussions with A. 
Brinkman, O. V. Dolgov, I. I. Mazin, U. Welp, A. Rydh, 
M. Iavarone, and G. Karapetrov. In Argonne this work 
was supported by the U.S. DOE, Office of Science, under 
contract # W-31-109-ENG-38. 



for the matrix elements used in Eqs. I|20(l . Then the tt- 
band F-function is given by 



V 2 z ( 2ttHx 
2 V *n 



(A2) 



Substituting this expression into the second self consis- 
tency equation, we represent it in the form 



-W21A1 + IU 2 A 2 



V2z ( 2nHx Y 
4vrT ^ $ J 



A, 



(A3) 

It has to be solved together with equations for F\ and 
the first self consistency equation. Using expansion with 
respect to eigenfunctions of the tr-band ^ n (x), this equa- 
tion reduces to the form of linear equation (|25b|l . in which 
the matrix Uij is given by the matrix elements 



U, 



dx^2i(x)^2j(x)g 



V 2z f 2irHx \' 
4^T \ $0 J 



Introducing the dimcnsionlcss oscillator wave functions 
tp n (x), ^n(x) = ^ n (x/xi)/ y/xl, we present these matrix 
elements in the dimensionless form 



APPENDIX A: LOCAL APPROXIMATION FOR 
THE 7T-BAND 

Let us consider equation for the F-function in the weak 
7r-band 



loF 2 



V 2z /2ttHx\- 
2 1 $ 



Fn 



2x VlF 2 



(Al) 



The typical scale of A2(a;) variation is imposed by the 
strong er-band. This scale is given by x\ = y -fi$>o/2TrH 
and, due to inequality 71 3> 72, it is much larger than the 
length scale X2 — ^72 < i ) o/2vriJ of the oscillator operator 
in the left side of the Eq. 1|A1J) . For relevant oj's the 
typical length scale of F2 variation is much larger than X2- 
This allows us to neglect the gradient term in Eq. i|Al|l . 
This approximation is equivalent to the approximation 



<^tp2i(x)ip2j{x)g [x /t 2 \ , 



where, again X/t2 = {T>2 Z /'D i z ) (H/tHi). In particular, 
ipo(x) = ti-" 1/4 exp(-i 2 /2). In the "high-field in 7r-band" 
regime, f2 < 1, one can use asymptotics g [x 2 /^] ~ 
2 ln(i ) — In £2 + 1e + 2 In 2 and obtain 

Ui,j(t 2 ) = fi,j + (ln< 2 - IE - 21n2)5„, m 

fi,j = -4 / dxip2i(x)^2j(x)lv-(x). 
Jo 

In particular, 

foo = — -7= / diexp(-£ 2 )ln(i) =7_E + 21n2 w 1.9635 

V 7T J 
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